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Exact large-iV/ results for the QCD Adler D-function and Deep Inelastic Scattering 
sum rules are used to resum to all orders the portion of QCD perturbative coefficients 
containing the highest power of 6=|(lliV-2iV/) , for SU(iV) QCD with Nf quark flavours. 
These terms correspond to renormalon singularities in the Borel plane and are expected 
asymptotically to dominate the coefficients to all orders in the 1/Nj expansion. Remark- 
ably, we note that this is already apparent in comparisons with the exact next-to-leading 
order (NLO) and next-to-NLO (NNLO) perturbative coefficients. The ultra-violet (UV) 
and infra-red (IR) renormalon singularities in the Borel transform are isolated and the 
Borel sum (principal value regulated for IR) performed. Resummed results are also ob- 
tained for the Minkowski quantities related to the -D-function, the e + e~ i?-ratio and the 
analogous r-lepton decay ratio, R T . The renormalization scheme dependence of these par- 
tial resummations is discussed and they are compared with the results from other groups 
[1-3] and with exact fixed order perturbation theory at NNLO. Prospects for improving 
the resummation by including more exact details of the Borel transform are considered. 
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1 Introduction 



In several recent papers [1-3] the possibility of resummation to all orders of the part of 
perturbative corrections contributed by QCD renormalons has been explored. 

The QCD perturbative corrections to some generic QCD Green's function or current 
correlator (to the Adler D-function of QCD vacuum polarization for instance) can be 
written: 

D = a + d x a 2 + d 2 a 3 H h d k a k+1 H , (1) 

where a = a s /ix is the renormalization group (RG) improved coupling; and the perturba- 
tive coefficients 4 can themselves be written as polynomials of degree k in the number 
of quark flavours, Nf] we shall assume massless quarks. 

4 = d fN k f + dt^N*- 1 + ■ ■ • + df . (2) 

The "iVy-expansion" coefficients, dj* r \ will consist of sums of multinomials in the ad- 
joint and fundamental Casimirs, C A =N, C F =(N 2 -1)/2N, of SU(iV) QCD; and will have 
the structure C\- r ~ s C F . The terms in this 11 N f -expansion" will correspond to Feynman 
diagrams with differing numbers of vacuum polarization loops. By explicit evaluation of 
diagrams with chains of such loops inserted, it has been possible to obtain the leading 4 
coefficient exactly to all orders for the Adler D-function [4-6] (and hence its Minkowski 
continuations, the e + e~ QCD _R-ratio and the r-decay ratio, R T ); the Gross Llewellyn- 
Smith (GLS) sum rule corrections ||; and heavy quark decay widths and pole masses 
A general procedure enabling d$ to be obtained from knowledge of the one-loop 
correction with a fictitious gluon mass has been developed |l|, 0. 

In a recent paper || we pointed out that the large order behaviour of perturbative 
coefficients is most transparently discussed in terms of an expansion of the perturbative 
coefficients in powers of 6=(11C^— 2Nf)/6, the first QCD beta-function coefficient: 

4 = 4^ + 4 fc - i V- l + --- + 4 0) . (3) 

This "6-expansion" is uniquely obtained by substituting Nf={^CA — 36) in equation (2). 
d\ ,= (— 1/3) fe 4 and so exact knowledge of the leading- Nf 4 to all orders implies exact 
knowledge of the 4 . 

In QCD one expects the large-order growth of perturbative coefficients to be driven by 
Borel plane singularities at z=Z£=2£/b, with £=±1, ±2, ±3, . . .. The singularities on the 
negative real axis are the so-called ultra-violet renormalons, UVg, and those on the positive 
real axis are the infra-red renormalons, I Re. These singularities result in the large-order 
behaviour of the coefficients 4 ~ h k k\. Indeed we showed in reference || that, given a set 
of renormalon singularities at the expected positions in the Borel plane, the leading terms 
in the 6-expansion, d^b k , should, if expanded in powers of Nf, asymptotically reproduce 
the 4f ^ coefficients of equation (2) up to 0(1/ k) accuracy. We conversely checked that 
the exact d^'b k results corresponded to a set of renormalon singularities at the expected 
positions. We shall demonstrate in section 2 of the present paper that, for the Adler 
.D-function and the GLS sum rule, the iVy-expansion coefficients obtained by expanding 
4 b and d^ b 2 are in good (10-20% level) agreement with those of the exact 0(a 3 ) next- 
to-next-to-leading order (NNLO) perturbative calculations for these quantities, so that 
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the anticipated asymptotic dominance of the leading-6 term is already apparent in low 
orders. 

Given the dominance of the leading-6 terms, an obvious proposal is to sum them to 
all orders. That is to split D into two components: 

D = D (L) + D {NL > (4) 

where 'L' and l NU superscripts refer to leading and non-leading terms in the fc-expansion. 

oo 

flW = £ 4% k a k+1 (5) 
fc=i 

and 

oo k— 1 

D^ = J2 ak+1 E d k } b e ■ (6) 

fe=l £=0 

The summation of terms can be achieved by using the Borel sum. The Borel integral can 
itself be split into two components and is well defined for the UVe singularities on the 
negative axis, which contribute poles to the Borel transform of D^ L \ The integral can 
be performed explicitly in terms of exponential integral functions and other elementary 
functions. The piece of the Borel integral for involving the I Re singularities on 
the positive real axis is formally divergent; but a principal value or other prescription 
can be used to go around the poles. The specification of this prescription is intimately 
linked to the procedure needed to combine the non-perturbative vacuum condensates in 
the operator product expansion (OPE) with the perturbation theory in order to arrive at 
a well-defined result for D ||. 

In recent papers by Neubert B and by Ball, Beneke and Braun [|I|, |2|] a summation 
of the leading-6 terms has also been considered. In these papers it has been motivated 
as a generalisation of the BLM scale fixing prescription (TO] and termed "naive non- 



abelianization" [II . The Neubert procedure uses weighted integrals over a running cou- 
pling. For the Euclidean Adler D-function this representation is equivalent to splitting 
the Borel integral into ultra-violet renormalon and infra-red renormalon singularities and 
principal value regulating the latter. When one continues to Minkowski space to obtain 
the e + e~ i?-ratio and the r-decay ratio, R T , there are several inequivalent ways to perform 
the continuation of the running coupling representation; and hence apparent additional 
non-perturbative ambiguities are claimed. In reference |2| the resummation is defined by 
using the principal value regulated Borel integral, as we shall do. They concentrate on R T 
and heavy quark pole masses. We agree with reference || that only a consideration of the 
singularities in the Borel integral provides a satisfactory way of combining perturbative 
effects with non-perturbative condensates, along the lines discussed in reference ||; and 
that the extra uncertainties claimed in reference || are spurious. 

Our intention in this paper is to focus on the Adler D-function, the e + e~ i?-ratio, R T 
and the GLS sum rule (the latter was not considered in references [1-3]). For all of these 
quantities there exist exact NNLO fixed order perturbative calculations and our interest is 
in comparing the leading-6 resummation with these exact fixed order results. The large-6 
results provide partial information about the Borel transform and the question is how 
this can best be utilised. We discuss the renormalization scheme (RS) dependence of 
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the split between and D^ NL ^ in equation (4), the relative contribution of being 
RS-dependent. This RS uncertainty needs to be kept in mind and is carefully discussed. 

The organisation of the paper is as follows. In section 2 we shall discuss the 6-expansion 
for the Adler D-function and the GLS sum rule and will consider the extent to which the 
dominance of the leading-6 term is RS-dependent. In section 3 the exact large-6 results 
are used to determine partially the Borel transforms for these quantities; and, having split 
them into UV and IR renormalon pieces, a resummation along the lines discussed above 
is performed. 

In section 4 the RS dependence of the resummed results is considered, the numerical 
results are presented and a comparison with the results of references || H and with the 
exact NNLO perturbative results is made. A discussion of the uncertainties and the 
prospects for improving the resummation by including more of the exact structure of the 
Borel transform is then undertaken. Section 5 contains overall conclusions. 



2 The 6-expansions for D and K 

We begin by defining the Adler D-function and the GLS sum rule, the quantities with 
which we shall be concerned in this paper. 

D(Q 2 ) is related to the vacuum polarization function, the correlator, U(Q 2 ), of two 
vector currents in the Euclidean region, 

- g^q 2 MQ 2 ) = igttH J d 4 xe^- x (o|T{j M (x) J„(o)}|o) , (7) 

by 

D(Q 2 ) = -\Q 2 ^{Q 2 ) > (8) 

with Q 2 =-q 2 > 0. 

In perturbation theory one has 
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D(Q 2 ) = d(R) £ Q} (l + \C F D) + \X Q f \ D . (9) 

Here d(R) is the dimension of the quark representation of the colour group, d(R)=N for 
SU(iV) QCD, and Qf denotes the quark charges, summed over the accessible flavours at a 

given energy. D denotes corrections of the "light-by-light" type which first enter at 0(a 3 ) 
and will be subleading in Nf] they are, it is to be hoped, small. D represents the QCD 
corrections to the zeroth order parton model result and has the form 

D = a + d x a 2 + d 2 a 3 H h d k a k+1 H . (10) 

We shall next consider two QCD deep inelastic scattering sum rules. The first is the 
polarized Bjorken sum rule (PBjSR): 



KpBj = f gl P ~ en (x, Q 2 )dx 

J 



9a 



9v 



l-\c F K) . (11) 
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Here K denotes the perturbative corrections to the zeroth order parton model sum rule, 

K = a + K x a 2 + K 2 c? + ■■■ + K k a k+l + ■■■ . (12) 
We can also consider the GLS sum rule, 

K GLS = \ f F?*" (x,<?)dx 
6 Jo 

i-^c F k + k^j . (13) 

The perturbative corrections, K, are the same as for the PBjSR; but there are additional 
corrections of "light-by-light" type, K, analogous to D of equation (9). These will similarly 
enter at 0(a 3 ) and be subleading in Nf and we shall assume once again that they are 
small. 

For both D and K the first two perturbative coefficients, d\, d 2 and Ki, K 2 , are known 
from the exact perturbative calculations [12-15]. We shall assume MS renormalization 
with renormalization scale fi=Q for the present; but will later discuss RS dependence 
more generally. 

Writing d% and d 2 expanded in Nf as in equation (2), the exact calculations give: 

* = (-^h + ^(f-T & K c - (14 > 

V162 27 s3 / f V 81 27 s3 9 s5 / i 

29 19. 10. \„ ^ 2 / 90445 2737 . 55 



+C F + — C 3 Cs )Nf + Ci[ Ca C 5 

V 96 6 S3 3^J f A \ 2592 108 S3 18 s5 

_ _ / 127 143^ 55 \ „ 2 ( 23 \ 
+C -^ ("Is - + T & ) + Ci (-32) ' (15) 

Here (3 and (5 are Riemann (^-functions. For later comparisons it will be useful to write 
these results numerically for SU(iV) QCD. 

.655iV + — J , (16) 

d 2 = .086iVj + Nf f-1.40iV - + (2.WN 2 - .661 - . (17) 

Expanding d\ and d 2 in b as in equation (3) gives 

dl = (^_ 2C3 ) 6+ |_^, (18) 

/151 19 \ )2 _ /31 5 A 5^ 
_ /29 19, in ,\, ^ 2 / 799 



32 2 au a V V 288 

^c^+^Uc.f.a. (l9) 



192 2 ' V r V 32 

As observed in reference [|J, the 6-expansion exhibits certain simplifications relative to 
that in Nf. In particular the £3, present in all orders of the iV/-expansion for d\, is present 
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only in the leading term in the 6-expansion. For d 2 the £5, present in all but the leading 
term in the iV/-expansion, is now present only in ■ In both cases the highest ^-function 
present cancels and is absent in the 'conformal' b — > limit [I5|; so cZj does not involve 
£ 3 and does not involve (5. Since the MS beta-function coefficients do not involve 
^-functions, this is presumably not an artefact of the particular RS chosen but may well 
be of more fundamental significance. It may ultimately be connected with the fact that 
in the 6 — ^ limit the z=2£/b UVe and IRi singularities in the Borel plane move off to 
infinity, leaving only instanton singularities. As discussed in reference || the (^-functions 
are intimately linked with the presence of renormalon singularities. It is amusing to notice 
that in the original NNLO result for d 2 [O, which was subsequently found to be in error 



D|, d^ does contain a non- vanishing — term for SU(3) QCD. If a fundamental result 




about the absence of ^-functions in the conformal limit could be established it would have 
enabled the incorrect result to have been dismissed at once. 

The corresponding results for the A^-expansion for the deep inelastic sum rules are: 

Ki = -- t Xf - {^-Ca-^-Cf) , (20) 
1 V648/ J \ 1296 2 ' 9*7 " ■ J V864 18 J 7 * 

/KAIT KK \ / 11/11 1 1 \ / 1 \ 

V 648 18 7 V 432 9 7 V327 v ; 
In numerical form, for SU(A r ) QCD, 

/ 438 \ 

K x = -.333A7+ M.48AT+— J , (22) 

K 2 = .177N} + Nf (-2.51N - + (a.53N 2 + .686 + . (23) 

Expanding in powers of b as in equation (3), 

Kl = " + {%-\ Cf ) • (24) 

V 72 / U44 2 W 9 S 7 V 288 6 S 7 

Notice that K\ does not contain £3; but a similar remark about the absence of (5 in K 2 ^ 
holds. 

We now wish to demonstrate that the leading term in the 6-expansion, when expanded 
in Nf, approximates the A^-expansion coefficients well, even in rather low orders. 
For d\ and d 2 we have 

d^b = .3456= -.115A^ / + .634A^, (26) 
4 2) 6 2 = .7766 2 = .086iVj - .9A8N f N + 2MN 2 . (27) 

The subleading, N, NfN and A^ 2 , coefficients approximate well in sign and magnitude 
those in the exact expressions in equations (16) and (17). The leading, Nf and Nj, 
coefficients of course agree exactly. 
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For K\ and K% we have 




b = -.333N f + 1.83N , 

1.596 2 = .177N 2 f - IMNfN + 5.37N 2 . 



(28) 
(29) 



The agreement with the exact N, NfN and N 2 coefficients in equations (22) and (23) is 
again rather good. 

We now turn to a consideration of the RS dependence of the Nf and b expansions. In 
variants of minimal subtraction, where the 1/e pole in dimensional regularization is sub- 
tracted along with an iVj-independent finite part, K, the QCD perturbative coefficients 
will have the form of polynomials in Nf as in equation (2). Modified minimal subtrac- 
tion (MS), corresponding to K = (In Ak-^e), with 7^=0.5722 . . ., Euler's constant, is most 
commonly employed. We can consider MS with renormalization scale fi=e u Q, where u 
is an iVy-independent number. The most general subtraction procedure which will result 
in perturbative coefficients polynomial in Nf, however, can be regarded as MS with scale 
[i=e u+v / b Q, where v is again iVj-independent. We shall refer to such renormalization 
schemes as 'regular' schemes. Of course the renormalization scheme is not specified by 
the scale and subtraction procedure alone but by higher order beta-function coefficients 
as well. Any variant of minimal subtraction with an iVf-independent renormalization 
scale will have v=0. Momentum space subtraction (MOM) based on the ggg vertex at a 
symmetric subtraction point fi 2 =Q 2 JTj| corresponds to u=2.56 and v=CAf{£,), where / 
is a cubic polynomial in the gauge parameter £. For the Landau gauge, £=0, v=— 2A9Ca- 
For other versions of MOM based on the qqg or ghost vertices, v will involve Ca and Cp- 

Let us denote the perturbative coefficients in the MS scheme with \i=Q (u—v—0) by 
dk] and those with general u and v by d' k . Then 



Changing v, one can make the d\ coefficient as large as one pleases and hence destroy 
the dominance of the leading-6 term noted above for the D-function and the sum rules in 
low orders; although the leading-6 term should still reproduce asymptotically the <ij|f 
coefficients to 0(1/ k) accuracy. 

For g?2 and higher coefficients the specification of the RS will involve higher beta- 
function coefficients as well as the scale and subtraction procedure. The RG-improved 
coupling, a(fj, 2 ), will evolve with renormalization scale according to the beta-function 
equation 



i 



(d^ + u )b+(dr +v) 

d\ + bu + v . 



(30) 




—ba 2 {l + ca + C2d 2 + 



• • + c k a k + •••)• 



(31) 



Here b and c are universal with [|TJJ 



b = hllC A 
o 

r 7c 2 A 



2N f ), 

11 C A C F 
8 b 



(32) 
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Integrating equation (31) with a suitable choice of boundary condition [pOfl , one obtains 
a transcendental equation for a: 



b In 



u 1 ca /" a , 

— = — hem h / ax 

A a 1 + ca Jo 



1 



x 2 B(x] 



+ 



cx) 



(33) 



The beta-function coefficients, 



where B(x) = (l + cx + C2X 2 + C3X 3 + • • • + CkX k + ■ ■ 
C2, C3, ... together with 6 In ^ label the RS. In a fixed order perturbative calculation one 
would truncate the beta-function. For the all-orders resummations of the next section, 
however, one requires an all-orders definition of the coupling. In the MS scheme the higher 



beta-function coefficients, d^ s , d^ s , . . ., presumably exhibit factorial growth, cj^ s ~ k\; 
and the 'a' coupling in the Borel integral would not be defined, since B(x) would itself 
need to be defined by a Borel integral or other summation. One therefore needs to use 
a finite scheme ]2]J where B(x) has a finite radius of convergence and can be summed. 
An extreme example is the so-called 't Hooft scheme [21] where 02=03=- ■ -=Ck=- ■ -=0, 
B(x) = l + cx. This results in the all-orders definition of the coupling, 



b In ^- = — h c In 

A a 1 + ca 



(34) 



In such a finite scheme, where C2, C3, ... are ^"/-independent, the 6-expansion of equa- 
tion (3) will contain an extra d k ^ /b term (for k > 1) and the dk will strictly no longer be 
polynomials in Nf. bdk, however, is a polynomial in Nf of degree k + 1. The leading-6 co- 
efficients, d^\ in regular schemes are independent of c 2 , c 3 , . . ., since these beta- function 
coefficients, Ck, are 0{1/Nf) relative to dk- 



3 Leading-6 Resummations 

Let us begin by recalling the definition of the Borel transform. If D has a series expansion 
in 'a' as in equation (1), we can write 

f°° 1 

D = / dze- z/a B[D](z) , (35) 
Jo 

where B[D](z) is the Borel transform of D, defined by (do=l) 

00 z m d 

B[D)(z) = £ . (36) 



m =o m - 

For the Euclidean quantities, D and K, defined earlier we shall deduce from the exact 
large- Nf results that, in the MS scheme with /z=e -5 / 6 Q, B[D](z) is of the form 

mm ) y> A Q (£) + A^)z + AijQz + M^ 2 + ■■■ 

~ B (£) + B x {t)z + B x {t)z + B 2 (£)z 2 + ■■■ 
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where Z£=2£/b. The two terms correspond to a summation over the ultra-violet renor- 
malons, UVe, and infra-red renormalons, IRi, respectively. Aq(£), A\(£), ae and B (£), 
B\(£), f3i will be obtained from the large- Nf results. The barred terms are sub-leading in 
Nf and remain unknown. The use of the so-called 'V-scheme' 0, MS with fi=e~ 5 / 6 Q, 
means that only the constant and O(z) terms in the numerator polynomials are leading 
in Nf. For a general MS scale, fi=e u Q, an overall factor e fez ("+ 5 / 6 ) should multiply the 
unbarred leading- Nf terms in the numerator. The presence of this exponential factor, 
when it is expanded in powers of z, can mask the presence of the UV and IR renormalons 
in low orders of perturbation theory. 

Notice that, whilst the residue at each renormalon singularity is only known to leading 
order in Nf, the A (£) and B (£) constant terms in the numerator polynomials are known 
exactly; indeed J2eLi(Ao(£) + B (£))—l, as is required to reproduce the unit coefficient of 
the 0(a) term in D in equation (1). We now turn to the explicit determination of the 
coefficients and exponents for D and K. 

For D the leading- Nf terms in the QED Gell-Mann-Low function are generated by || 



where 




=1 



32 ~ (-l) k k , , 

PW = 3(iT^)£(kW- (38) 

In the V-scheme, MS with n=e~ 5 ^ 6 Q, one then has 

#>=(~) n 2^. (39) 

It is then straightforward to deduce that the coefficients and exponents in equation (37) 
for B[D](z) are 



3 £ 2 (£+l) 2 (£ + 2) 2 ' iW 3£ 2 (£+l) 2 (£ + 2) 2 
£ = 1,2,3,... 

B (1)=0, B (2) = l, B (£) = -A (-£) £>3 
B 1 (1)=0, Si(2) = 0, B x (£) = -Aii-i) £>3 

a e = 2 £=1,2,3,..., p 2 = l, fa = 2 £>3. (40) 

So IR\ is absent, as required from the absence of a dimension two condensate in the 
OPE 0, ||. IR 2 is a single pole. All the other singularities are double poles. Not only 
are the coefficients for the UVi and IRi singularities related by the curious symmetry 
Bo,i{£)=— Aqi(— £); but the form of A (£) means that there is an additional relation, 
Ao{£)=— Bq(£+2), so that the constant term in the numerator polynomial for UVe exactly 
cancels that for IRi +2 - This ensures that 

oo 

Y / (Ao(£) + B (l)) = B (2) = l, 
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n 




UVi 


UV 


IR 





1 


.8148148 


.7198242 


.28018 


1 


-.4874471 


-.6296296 


-.5921448 


.10470 


2 


.8938293 


.8518519 


.8258924 


.06794 


3 


-1.525257 


-1.611111 


-1.586113 


.06086 


4 


3.927235 


3.888889 


3.858335 


.06890 


5 


-11.24973 


-11.38889 


-11.34378 


.09405 


6 


39.23893 


39.16667 


39.08871 


.15022 


7 


-154.1541 


-154.5833 


-154.4291 


.27499 


8 


688.5574 


688.3333 


687.9894 


.56801 


9 


-3410.339 


-3412.500 


-3411.647 


1.3079 


10 


18638.50 


18637.50 


18635.17 


3.3243 



Table 1: Leading-6 coefficients, for the Adler D-function, D, compared with the 
contribution of the first UV renormalon, { UV\\ (equation (41)). The V-scheme, MS with 
[i=e~ 5 / 6 Q, is assumed. l UV and l IR denote the separate sums over the UVi and IRi 
singularities. 



which, as noted earlier, is required to reproduce the unit coefficient of the 0(a) term in 
the perturbative expansion. The precise origin of these relations between UV and IR 
renormalons remains unclear and deserves further study. They have also been noted and 



discussed in reference 22 



For the D-function the singularity nearest the origin is UV\ and from the Ao(l), -Ai(l) 
in equation (40) this should correspond to 

, (n)| 12n + 22 / l\ n 

drluvi = s= — n\ [--) . (41) 



27 V 2 

In Table 1 we compare the exact leading-6, coefficients with the contribution 
from UV\ of equation (41). UV and IR denote the separate sums over the UVi and IRi 
singularities. The MS scheme with fi=e~ 5 ^ 6 Q (V-scheme) is assumed. With this choice 
of scheme UV\ dominates even in low orders and the alternating factorial behaviour is 
apparent. 

We now consider the coefficients and exponents in equation (37) for B[K\(z). The 
generating function for the leading-6 coefficient in the V-scheme is M 



K (n) 

n 

This results in 



1 /l\ n d n (3 + x) 



3 \2J dx"(l-x 2 )(l-^; 



(42) 

x=0 



4 18 5 



The terms correspond to UV±, UV2, IRi, IR2 respectively. Each numerator and 
exponent will contain in addition 0(1/Nf) corrections corresponding to the barred terms 
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n 


n 


UV l + IRi 


UV 


IR 





l 


1.333333 


.3888889 


.6111111 


1 


.1666667 


.2222222 


-.2083333 


.3750000 


2 


.6250000 


.6666667 


.2152778 


.4097222 


3 


.3125000 


.3333333 


-.3281250 


.6406250 


4 


1.968750 


2.000000 


.6614583 


1.307292 


5 


1.640625 


1.666667 


-1.660156 


3.300781 


6 


14.94141 


15.00000 


4.990234 


9.951172 


7 


17.43164 


17.50000 


-17.48291 


34.91455 


8 


209.7949 


210.0000 


69.96582 


139.8291 


9 


314.6924 


315.0000 


-314.9231 


629.6155 


10 


4723.846 


4725.000 


1574.808 


3149.039 



Table 2: As for Table 1 but for the Deep Inelastic Scattering sum rules, K, 
t UVx+IRi denotes the contribution of the singularities nearest the origin, the first two 
terms of equation (44). 



in equation (37). The constant terms in the numerators sum to 1, again ensuring a unit 
0(a) coefficient in K. It would be interesting to try to understand the fact that only the 
first two UV and IR renormalons are leading in Nf in the context of the OPE for the 
deep inelastic sum rules, a topic discussed in reference |23]. 
is then given by (in the V-scheme) 



9 



9 



-ni - 



-m 



(44) 



Table 2 shows that K^y is dominated, even in low orders, by the combined UV% + IR\ 
contributions of the two singularities nearest the origin, the first two terms of equation 
(44). 

We now wish to use the Borel integrals to perform the leading-6 resummation defined 
in equations (5,6). For the Adler D-function we have 



D {L \a) 



A {£)+A 1 {£)z 




(45) 



The coefficients Aq, A\, Bo, B\ are summarised in equation (40). We assume that the 
resummation has been performed in a finite RS corresponding to MS subtraction with 
[i=e u+v / b Q. 'a', the coupling, is then defined by the integrated beta-function equation of 
equation (33). One then finds that the exponent in equation (45) is 



F(a) = bin 



Q 



A 



MS 



5 ca 
—b — cm 

6 1 + ca 



+ v — / dx 



+ 



x 2 B(x) x 2 (l + cx) 



(46) 
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The RS dependence of D^ L > reflects the fact that only a subset of the perturbation series 
has been resummed, hence violating the exact RS-invariance which would apply to the 
full series. We shall return to this RS dependence in a moment. 

The first, UV renormalon, term in equation (45) is a completely well-defined inte- 
gral. It may be performed in terms of the exponential integral function (with negative 
argument), 

/OO p — t 

dt — . (47) 
-X t 

The first term yields 

oo 

D {L \a)\ uv = ]>>{e F(a) ^Ei(-F(a)^) [F(a) Zi (A(£) - z,A x {t)) - z t A x {i)\ 

+(A (£) - z t A x {l))} . (48) 

To evaluate the second, IR renormalon, term we shall use a principal value prescription; 
correspondingly we need to define Ei(x) with a positive argument as a principal value. 
We find 

D {L \a)\ m = e- F ^ Z2 z 2 B (2)Ei(F(a)z 2 ) 

oo 

+ E^{e" F(a)2< Ei(F(a)^) [F(a)z e (B (£) + z e B 1 (£)) - zeB^i)] 

£=3 

-(B (£) + z e B ± (£))}. (49) 

Finally 

D^(a) = D^(a)\uv + D^(a)\ IR . (50) 

For the Deep Inelastic Scattering sum rules we will have, analogously, using equation 

(43), 



~ poo 

K {L \a) = / d^e- F(a) 
Jo 



4 J_ 

9 18 



(i + t) + 
dze -^(_i___A_). (51) 

These integrals may be expressed once again in terms of Ei(x): 

K (L) (a)\ uv = Ll e F ^ 1 z 1 Ei(-F(a)z 1 ) + ±- e F(a >>z 2 Ei(-F(a)z 2 )} (52) 

y to 

and 

K (L \a)\ IR = ^e- F ^ Zl Ei(F(a) Zl ) - ^e~ F ^ z 2 Ei(F(a)z 2 )j . (53) 

Similarly 

K^ L \a) = K^ L \a)\ uv + K^ L \a)\ IR . (54) 

Before we numerically evaluate these results and comment further on RS dependence, 
we shall derive the analogous resummations for the Minkowski continuations of the IP- 
function, the e + e~ annihilation I?-ratio and the analogous quantity in r-decay, R T . The 
I?-ratio is related to IP by a dispersion relation, 
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Here s is the physical timelike Minkowski squared momentum transfer. A perturbative 
result for R of the form of equation (9) can be written down involving a quantity R 
with perturbative coefficients r^. The are directly related to the dk via the dispersion 
relation (55): r±=di, T2=d<i — n 2 b 2 /12. The 7r 2 term arises due to analytical continuation. 
In the Borel plane one finds (to leading order in Nf) 

B[R](z) = ^BP|(,) • (56) 
The leading-6 resummation is then obtained from equation (45) simply by adding an extra 



, factor in the integrand. 



sin(7rfrz/2) 

iauui in unc 111uc51a.il 

(a)z Bm(irbz/2)~MQ + MQz 



[ ) ~ Jo nbz/2 h + i 



z_\2 



22 ' e=3 



Writing the 'sin' as a sum of complex exponentials and using partial fractions, the UV 
integrals can be explicitly performed in terms of the generalised exponential integral 
functions Ei(n, w), with complex argument w, defined for Rew > by 



00 n—wt 



One also needs 



Ei(n,w) = I &t e —. (58) 



BY„i» sin(7rte/2) / irb \ , . 

dze~ F{a)z — i !—L = arctari — — . (59) 



z \ 2F ( a ) 



One finds 



xvn, \ 1 2 /8C2 H\ / vrfe 



r> OO 

+ ~ h ^\ + (4>(4 - A 1 (£)^)0 + (2,£) \ , (60) 



where 

+ (p,g) = e F ^(-iyim[Ei(p,F + z q )\ (61) 

with F±=F(a) ± if*. 

To evaluate the principal value of the IR contribution in equation (57) one needs to 
continue Ei(n, w), defined by equation (58) for Rew > 0, to Rew < 0. With the stan- 
dard continuation one then arrives at a function analytic everywhere in the cut complex 
u>-plane, except at w=0; and with a branch cut running along the negative real axis. 
Explicitly |32| 



Ei(n, w) 



{n- 1)1 



n—l -y 



■Into - 7 B + J2 — 



m=i m . 



„ (m — n + ljm! 

m=0 v ' 
m^n— 1 
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with 7^=0.5722 . . ., Euler's constant. The In w term in equation (62) means that Ei(n, w) 
is not a real function. For instance, for negative real w one has Ei(l, — x ±ie)=Ei(x) =F iir, 
where Ei(x) is the principal value of equation (47) used to define the IR renormalon 
contribution for the Euclidean quantities. 

In order to evaluate the IR renormalon contribution correctly, one in fact needs to 
continue Ei(n, w) as a real function, so that for Hew < one makes the replacement 
law — > lnu> + in sign(Im w). Correspondingly, one should define the IR analogue of 
equation (61), 



0_(p, g ) = e~ F ^(-l)nm[Ei(p, -F +Zq )] - 6 ^ M^)'" 1 ] 



(63) 



where Ei(p, —F + z q ) is defined by equation (62). The principal value of the IR renormalon 
contribution is then given by 



a)\m 



nb { 3 3 



arctan 



Tib 



2F(a) 



+ 



2B (2) 
Tib 



1,2) 



2 00 r i 

+^E{ 5 oW<MM) + (B (e) + B 1 (£)z e )4>^2,e)j. 



Then 



& L \a) = R {L \a)\ uv + R {L \a)\ IR 



(64) 



(65) 



The r-decay analogue of the i?-ratio, R T , can be defined in terms of the .R-ratio by 
the integral representation 



R T 



m t ds 

— (l-,/M T 2 ) 2 (l + 2 S /M T 2 ) J R( S ) 



d{R)(\V t 



ltd 



IK. 



1 + -C F R T 



Here R(s) denotes R(s) with the J2fQ} replaced by \V u d\ 2 + \ V US \ 2 
KM mixing matrix elements. R T has the form 



(66) 

1, where the V^'s are 



Rt 



a + r[a 2 + rja 3 + 



+ r T k a k+1 + 



(67) 



It is then straightforward to show that 



B[Rr)(z) 



sin(nbz/2) 
nbz/2 



bz_- 

2 / 



(1 



bz 



+ 



1 



8 • 



B[D]{z) . 



(68) 



Proceeding in a manner analogous to that for the i?-ratio, we find 

Tib \ 



Rf\a) 



uv 



2_ n%2 

Tib \ 3 
4 



11 \ 

— arctan 



2F(a) 



+— J2UomG(£) + H(e)) - ^A(£)G(£))0 + (1, 
nb e=1 

+F(£)(AoW-zMiW)0+(2,£)] 



(69) 
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and 



iVwIifl = — (- — I arctan I 




14 64, , , /23 32, 

H In 2 - 8C 3 + bzxi ln2 

3 3 SJ 1 V 3 3 

■^ + 641n2 + te 3 (^-^ln2))0_(l,3) 

128, , /2035 16, 
In 2 + bz 4 { + — In 2 



7776 81 



<M2,3) + 



±1 


f 247 


5bz±\ 


Tib 1 


^648 ~ 


162 J 


±1 


f 13 - 




nb 1 


I 432 


1728 ) 




00 



where 



+^Y^[(B (£)(G(-£) + H(-£))+z e B 1 (£)G(-£))cj ) 4l,£) 
+H(-£)(B {£)+ztB 1 {£))<f>_(2,£)) , (70) 



Then, as before, 



(£ + l) 2 (£ + 3) 2 (£ + 4) 2 ' w (£ + l)(£ + 3)(£ + 4) ' 

i?^(a) = flj. £ )(a)|uv + ^ L) (a)|/ii • (71) 



Before we proceed to discuss RS dependence further and to evaluate numerically these 
resummed expressions, we would like to make some remarks. The first concerns the ease 
of evaluation of both the Euclidean resummed expressions, equations (48), (49) and (52), 
(53), and those for the Minkowski quantities, equations (60), (64) and (69), (70). Even 
though these expressions contain infinite summations over the contributions for the UVe 
and IRi singularities, successive terms in the sums are strongly damped, with the result 
that, in order to obtain the three significant figure accuracy of the resummed results to be 
tabulated in the next section, it is only necessary to retain terms up to and including £=7 
in each sum. The resummations can then be straightforwardly and rapidly evaluated. 

The second remark concerns the connection between these explicit expressions for 
the principal value of the Borel sum and the inequivalent continuations of the running 
coupling representation to the Minkowski region in reference |J. It is straightforward to 
show that procedure '1' of reference for R and R r corresponds exactly to evaluating 
equations (60), (64) and (69), (70) using q) defined by equation (63) with the second 
term omitted, i.e. using the standard continuation of Ei(n, w) defined in equation (62). 
This does not produce the principal value of the Borel sum. Worse still, the Borel sum 
contains pieces involving single IR renormalon poles together with a sin ^ factor, which 
are well-defined and finite due to the compensating zero contained in the 'sin'. These 
contributions, which do not require regulation, are evaluated incorrectly with the standard 
continuation. Procedure '2' of reference corresponds to evaluating the Borel sum, 
incorrectly omitting the second term in equation (63) for some of the IR singularities 
and, correctly, retaining it for others. In our view the inequivalent continuations of the 
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running coupling representation of reference || to the Minkowski region correspond to 
various ways of wrongly evaluating the Borel sum. We see no reason to believe that 
these discrepancies have a physical relevance, or that they reflect inadequacies in the 
definition of the OPE in the Minkowski region, as suggested in reference ||. We agree 
with reference @ that, with our present state of knowledge, the regulated Borel sum 
provides a satisfactory framework for combining IR renormalon ambiguities with the 
vacuum condensate ambiguities in the OPE. 



4 RS-dependence of the resummed results 



Armed with these resummed expressions, we now return to the question of the RS depen- 
dence of D^ L \a) before presenting the numerical results. 
For a generic quantity we can write 



D = £> (i) ( 



D^ NL \a) 



(72) 



'£)' on the left of the equation denotes the full Borel sum, i.e. equation (37) with barred 
and unbarred terms included and the IR singularities principal value regulated; and we 
assume that this exists. Crucially, the full Borel sum is RS independent and so will not 
depend on 'a'. The L and NL components do depend on 'a', however. We can consider 'a' 
varying between a=0 and a=+oo, labelling possible RS's. From equation (46), as a — > 
+oo, resulting from the — cln term, and we have assumed c > 0, which is 



so F(a) 



true for Nf < 8 in SU(3) QCD. One then has D^ L \0)=0. Correspondingly, from equation 
(72), D^ NL \0)=D and so, as a — > 0, the NL component contributes the whole resummed 
D. As 'a' increases the — cln^|^ term in equation (46) decreases and as a — > oo it 
vanishes, resulting in a finite limit F(oo): 



F(oo) = bin 



Q 



A 



MS 



6 



dx 



1 



x 2 B(x) x 2 (l + cx) 



(73) 



We assume that B(x) is such that the integral exists. Thus D^ L \a) increases from 
^^(O)^ to a finite maximum value, D^(oo), as 'a' increases; correspondingly, D^ NL \a), 
which provides the whole resummed D at a—0, decreases as 'a' increases. 

This RS dependence of D^ L \a) is clearly problematic. It is monotonic and hence there 
is no basis for choosing a particular scheme. There is also a dependence on the particular 
finite scheme, characterised by the choice of B(x), and on the parameter v. The maximum 
value, D( L \oo), does perhaps minimize the relative contribution of the unknown D^ NL ^ 
component but there is no guarantee that D^ NL \oo) is positive; and it is entirely possible 
that D^ L \oo) overestimates D. 

We shall choose v—0, corresponding to a variant of minimal subtraction, a choice 
which one can motivate by the observed dominance of the leading-6 term in MS noted in 
section 2. For reasons of simplicity we shall choose the 't Hooft scheme, corresponding to 
B(x)—1 + cx. With these choices one has 



F(oo) = bin 



Q 



A 



(74) 



MS 
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and we shall use this in the resummations. It corresponds simply to taking 'a' as the 
one-loop coupling in the MS scheme with n=e~ 5 ^ 6 Q (the V-scheme): 

where Ay=e 5 / 6 A^jg. This is in fact the same choice for a as in references [1-3], where it is 
motivated by noting that using the one-loop form for a(/x 2 ) makes the leading-6 summation 
/^-independent. We stress once again that in our view its significance is that it maximizes 
D^ L \a) for a given choice of finite scheme, B(x), and parameter v. In Figure 1 we show 
D^ L \a) plotted versus 'a' ('t Hooft scheme with v—0). In the figures we have plotted 
versus a = 1 — e~ a , so that the full RS variation can be fitted in a unit interval in a. We 
have taken Q=M z =91GeV and Aj^g(iVy-=5) = lllMeV. The solid curve gives the overall 
L)( L \a), split into D^ L \a)\uv (dashed) and D^ l \o)\ir (dashed-dot) contributions. Similar 
curves for K^ L \a) with <5 2 =2.5GeV 2 and A(^g(A^/=3)=201MeV are given in Figure 2; and 
for w L \a) with Q=91GeV in Figure 3. The corresponding curve for w^\a) is given in 
Figure 4, Q=M T =1.78GeV and A M g(A^ / =3)=201MeV. The qualitative behaviour is as we 
described earlier. The relative sizes of the UV and IR contributions reflect the disposition 
of the UV and IR singularities described above for the different quantities. 
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0.015 L 
0.010 
0.005 _ 



Figure 1: The leading-6 resummation, D^ L \a), plotted versus a=l — e~ a ('t Hooft scheme 
v=Q), Q=M z =91GeV and Aj^g(iV/=5) is as in the text. The solid curve is the overall 
result split into D^ L \a)\uv (dashed) and D^ L \a)\ IR (dashed-dot) contributions. 
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Figure 2: As for Figure 1 but for K^^a) with Q 2 =2.5GeV 2 . 
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a 

Figure 3: As for Figure 1 but for R^(a) with Q=M z =91GeV. 

The a — > oo limits obtained using F(oo) as in equation (74) are tabulated in Table 3. 
Specifically, the 'Resummed' column contains l + D^ L >(oo), 1 + R( L \oo), l + R[ L \oo) and 
1 — i^ L )(oo) for each observable and we have added some extra energies, Qo=20GeV and 
Q?=2.5GeV 2 . The values of A^g for N f =5, 3 are as noted above. In each case we have 
taken care to include sufficient terms in the summation over UVe and IRe singularities to 
guarantee accuracy to the quoted number of significant figures. 

As noted, using F(oo) is equivalent to the one-loop definition of the coupling used in 
references [1-3]; and the same value of Aj^g(Afy=3) has been used. Our resummed results 
for D(m^), R(w%) agree with the principal value of the Borel sum for these quantities 
quoted in reference ||; and the resummed result for R T agrees with the results quoted 
in references [@, [|. As discussed at the end of section 3 we can reproduce the results of 
procedures '1' and '2' for R and R T in reference || by incorrectly omitting the second 
term in equation (63) for some of the IR renormalon singularities. 
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Observable 


Energy 


Resummed 


FOPT 


Expt 


D 


m T 


1.151 


1.087 


— 




Qo 


1.055 


1.045 


— 




m z 


1.042 


1.035 




R 




1.105 


1.080 






Qo 


1.053 


1.044 






m z 


1.041 


1.035 


1.040 ± 0.004 


R T 


m T 


1.228 


1.115 


1.183 ±0.010 


K 


Qi 


.784 


.889 


.768 ±0.09 



Table 3: Comparison of resummed results ('Resummed') of section 3 using F(oo), equation 
(74) (see text). Q 2 =(20GeV) 2 , Q 2 =2.5GeV 2 . 'FOPT' gives the exact NNLO perturbative 
results with fi=Q MS scheme and the values noted in the text. 'Expt' gives the 
experimentally deduced values for some of these quantities [25-27]. 



There is clearly an ambiguity associated with the IR renormalons or, correspondingly, 
with vacuum condensates in the OPE. For a leading IRg singularity in which there is a 
single pole one would expect an ambiguity in the principal value ~ B e~ Ze ^ a , where B is the 
residue of the renormalon. For D the leading IR singularity is /i?2- For the Minkowski 
quantities, R and R T , the sm 7 [^'^ 2 ' > factor apparently removes the single pole at z=z 2 
but there is presumably still a branch point singularity at z=z% beyond the leading-iVy 
approximation ||. The determination of the residue B would require a resummation of 
the numerator polynomial and so it is only known to leading- Afy (for D Bq(2)=1). Taking 
B=l and putting Nf=3, 5 values for b and 'a' values corresponding to the energies and Aj^g 
considered in Table 3 yields an ambiguity <10~ 4 for D, R, R T , so the significant figures 
quoted in the resummed result do not change. For K, however, the leading singularity is 
IRi and one would estimate the IR ambiguity ~ 10~ 2 for Nf=3 and Q 2 =2.5GeV 2 , which 
is clearly significant. 

The column labelled 'FOPT' gives the fixed order perturbation theory results ob- 
tained at NNLO (up to 0(a 3 )) using the exact perturbative coefficients in the MS scheme 
with fi=Q [12-15]. The coupling 'a' is defined using the NNLO truncated beta-function, 
B(x) = l + cx + c^ s x 2 , in equation (33), with the values of A^jg as above. 

The column labelled 'Expt' gives the values determined from experimental data for 
(1+R) [|25|, {1+R T ) @ and (1-K) |27|. The results of adjusting A^g to fit the FOPT and 
resummed predictions to these experimental values are summarised in Table 4. Both fixed 
order perturbation theory and the leading-6 resummed results exhibit RS-dependence and 
it is not at all obvious which procedure gives the closest approximation to the all-orders 
sum. We will defer a more detailed discussion of this question, including a consideration 



of how the resummed results compare with use of the effective charge formalism |28], ^9 
until a future work. 

To conclude this section let us consider how the leading-6 resummation might be 
improved. An obvious improvement would be to include the full branch point structure 
of the renormalon singularities by incorporating the subleading in Nf, a>e and /3 e , pieces 
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Observable 


N f 


Aj^g/MeV fitted to experimental data 






Resummed 


NNLO FOPT 


R 


5 


98t| 


287±?1| 


R T 


3 


159l? 


386t| 


K 


3 


37511 


495lil 



Table 4: Values of A^g adjusted to fit the predictions of NNLO FOPT and the resummed 
results to the experimental data for (1+R) PJ, (1+R T ) P§ and (1-K) [E7|]. 



of the exponents in equation (37) into the resummations. One expects 

a e = -cz e + it , 

P e = cz e + if, . (76) 

The first term can be deduced from RG considerations but the 7^ and 7^ are the one-loop 
anomalous dimensions of the relevant operators f9|, [30[ . For D it is known that IR2 has a 
corresponding OPE operator with vanishing one-loop anomalous dimension, 72=0 [^TJ, so 



(3 2 =cz2 is known. To the best of our knowledge the remaining 7^ and 7^'s are not known. 
One could nonetheless include the first RG-predictable terms in equation (76) in the 
resummations and see by how much the results change. The problematic RS-dependence 
of D^ L \a) would be qualitatively unchanged, however. 

Further improvement could be achieved by including some of the subleading in Nf, 
barred, coefficients in the numerator polynomials in equation (37). A complete fixed 
order perturbative calculation for D up to 0(a n ) would enable the series coefficients 
of B[D](z) up to 0(z n ) to be determined exactly; but to obtain the coefficients of the 
numerator polynomials for each singularity up to 0(z n ), even given knowledge of the full 
branch point exponents discussed above, would still be very difficult. The improvement 
of perturbation series by developing a representation of the form of equation (37) with 



truncated numerator polynomials was suggested in reference |3T . 

One further caveat concerns the structure of the UVg singularities. It has been sug- 
gested in reference |$0| that diagrams with more than one renormalon chain will modify 
the form of the ultraviolet renormalon singularities ||; this would lead to additional UV 
terms in equation (37) with O(z) or higher, subleading in Nf, numerator coefficients but 
with increasing leading- Nf exponents ae. The presence of such terms would destroy the 
asymptotic dominance of the leading-6 coefficient to all-orders in Nf 0, which seemed 
to be already evident in the comparisons with the exact NLO and NNLO coefficients 



23 



discussed in section 2. The motivation for leading-6 resummation would hence disappear 
if this UV structure is correct. Further clarification of this point is obviously required. 

5 Conclusions 

In this paper we have investigated the possibility of resummation to all orders of the 
leading term in the '6-expansion' of QCD perturbative coefficients in equation (3). This 
expansion was introduced and motivated in reference by a consideration of renormalon 
singularities in the Borel plane; and, if such singularities are present, then the d^b k term 
should, when expanded in Nf, reproduce the iV/-expansion coefficients of equation (2) 
to all orders in Nf with asymptotic accuracy 0(l/k). We checked explicitly in section 2 
that for the QCD Adler D-function (D) and Deep Inelastic Scattering sum rules (K) this 
asymptotic dominance of the leading-6 terms was already evident in comparisons with 
the exact NLO and NNLO perturbative coefficients for those quantities. The interesting 
absence of ^-functions from the d% coefficient which survives in the b — > limit was also 
noted. The RS dependence of the leading-6 coefficient and the need to give an all-orders 
definition of the coupling in order to perform resummations was discussed. 

In section 3 we used exact large-iVy results [4-6] to obtain partial information about 
the Borel transforms B[D](z) and B[K](z). Ultra-violet and infra-red renormalon singu- 
larities are present and we obtained the constant coefficients in the numerator polynomials 
exactly; and the exponents, single and double poles, to leading- iVj (equations (40), (43)). 
We showed that in the V-scheme (MS with fi=e~ 5 ^ 6 Q) the leading-6 coefficients, 
and KM, are dominated, even in low orders, by the renormalon singularities nearest the 
origin, respectively UV\ and UV\+IR\ combined (see Tables 1 & 2). 

For each quantity we split the leading-6 Borel sum into UV and IR poles. The first 
contribution could be evaluated exactly in terms of the exponential integral function 
Ei(x) (equation (47)) and elementary functions; and the second, IR, contribution could 
be obtained as a principal value, in terms of a principal value of Ei(x). We showed how 
to modify the Borel transform for the Minkowski continuations of D, the e + e~ i?-ratio 
(R) and the i?-ratio for r-decay (R T ), and a similar resummation was performed for these 
Minkowski quantities in terms of a generalised Ei(n, w) function (equations (58) and (62)). 

In this way we obtained the component of the split defined in equation (4) for 
the above quantities. Unfortunately the result obtained by summing the leading-6 terms 
is RS-dependent, D^ L \a). In section 4 we showed that one can maximize D^ L \a) and 
hence perhaps minimize D^ NL \a) for any given choice of finite scheme and subtraction 
procedure. Maximizing D^ L \a) whilst choosing the 't Hooft scheme and a variant of 
minimal subtraction (v=0) was equivalent to using the one-loop coupling in the V-scheme, 
the choice also made in references [1-3]. 

We compared our resummed results with the principal values of the Borel sum for 
D(m%) and R{m*) quoted in reference || and with that for R T in references [0, |I| and 
found agreement. They are tabulated in Table 3. The procedures '1' and '2' for continuing 
the running coupling representation of reference |J, for R and R T , to the Minkowski 
region were shown to correspond to using different continuations of the Ei(n, w) function 
from Hew > to Kew < 0. Only a continuation of Ei(n, w) as a real function enables 
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one to evaluate correctly the well-defined pieces of the Borel sum involving single IR 
renormalon poles with a sin ^rp factor. The inequivalent procedures of reference || are 
seen to correspond to various ways of wrongly evaluating the Borel sum and are therefore 
spurious. We stress that the Borel sum with IR singularities identified and principal 
value regulated provides a unique result and, in our view, a firm foundation for combining 
vacuum condensates in the OPE with IR renormalons to achieve a well-defined overall 
result (see reference 0). There is, of course, ambiguity due to the IR poles and we 
estimated this to be <10~ 4 for the D, R, R T resummed results in Table 3, so the displayed 
significant figures should be valid; but much larger, ~ 10~ 2 , for K due to the presence of 
an IRi singularity. 

We also compared with fixed order perturbation theory up to NNLO for these quan- 
tities. Since both D^ L \a) and the NNLO fixed order result suffer from RS dependence it 
is unclear which is the more reliable and further investigation of this question is required. 

We finally considered how the leading-6 resummation might be improved by including 
more exact information about the Borel transform. 
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